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ABSTBACT 

iiaum  freguency  for  propagation  of  sound  in  shallow 
sound  channels  was  studied  using  two  acoustic  transnission 
loss  oodels.  The  split-step  Parabolic  Eguation  model  (a 
full-wave  model)  and  the  Fast  Asymptotic  Coherent 
Transmission  loss  model  ^‘-/version  9H  (a  ray-tracing  model)  ?  - 
were  tested  against  experimental  data  collected  by  Dosso  and 
Chapman  in  the  northeast  Pacific  Ocean.  The  models  were 
found  to  be  valid  predictors  of  optimum  freguency  for  the 
shallow  sound  channel  observed  by  Dosso  and  Chapman.  Both 
models  were  then  used  to  predict  optimum  freguency  for  two 
sound  velocity  profiles  obtained  in  a  high-latitude  deep 
ocean  basin  under  summer  conditions,  exhibiting  shallow 
sound  channels.  As  expected,  the  split-step  Parabolic 
Equation  (PE)  model  adequately  predicted  optimum  frequencies 
for  these  cases.  The  Fast  Asymptotic  Coherent  Transmission 
loss  model,  version  9H  (FACT  9H)  model  did  not  produce 
reasonable  results  for  optimum  frequencies. 
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III.  THE  DCSSO  AND  CHAP HAN  EXPEBIHENT 


A.  DATA  COLLECTION 

Dcsso  and  Chapman  [Hef.  1]  conducted  an  experiment  in 
sound  propagation  within  and  below  the  shallow  sound 
channel.  During  summer  conditions  in  the  northeast  Pacific 
Ccean  off  the  west  coast  of  Canada,  CFAV  ENDEAVOUR  took  up  a 
station  and  monitored  data  collection  equipment.  A  single 
calibrated  hydrophone  was  suspended  below  the  vessel  at  a 
depth  of  101  meters,  approximately  the  shallow  sound  channel 
axis.  An  array  of  hydrophones  was  suspended  at  417  meters 
near  the  deep  sound  channel  axis.  A  second  vessel,  CVAF 
PARIZEAU,  opened  range  from  CFAV  ENDEAVOUR,  dropping  explo¬ 
sive  charges  every  1.8  km  (approximately  one  nautical  mile). 
The  charges  were  set  to  detonate  at  98  meters,  near  the 
shallow  sound  channel  axis. 

CFAV  PARIZEAU  also  measured  sound  velocity  profiles  at 
16-kilometer  intervals  over  the  track  (Figures  3.1  through 
3.3).  later  analysis  revealed  that  the  sound  channel  deep¬ 
ened  over  the  portion  of  the  track  between  eight  and  twenty 
kilometers,  so  that  the  sources  were  not  located  within  the 
SSC  and  sound  trapping  was  greatly  reduced  for  sources  deto¬ 
nated  between  those  ranges.  Section  D  below  discusses  steps 
that  Dosso  and  Chapman  took  to  account  for  this  deepening  of 
the  SSC. 

The  received  acoustic  energy  was  measured  as  1/3  octave 
band  averages  for  each  shot  for  both  the  shallow  and  deep 
hydrophones.  Dosso  and  Chapman  were  able  to  distinguish 
bottom-interacting  propagation  paths  by  the  different 
arrival  times  for  shcts  in  the  first  35  kilometers  of  the 
track,  and  these  arrivals  were  subtracted  from  the  resulting 
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Figure  2.2  Sound  Telocity  Profile 
Shoving  Shallow  Sound  Channel. 
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Figure  2.1  Typical  Hid-latitude  Sound  Velocity  Profile. 
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and  below  by  smooth  interfaces  between  layers  with  no 
absorption  or  leakace  should  produce  losses  at  about  the 
rate  due  to_  cylindrical  spreading  alone.  The  minimum 
frequency  which  will  be  trapped  well  in  a  channel  is  the 
cutoff  frequency.  At  frequencies  near  1000  Hz,  absorption 
mechanisms  become  important,  and  scattering  losses  also  come 
into  play.  Therefore  it  should  be  possible  to  find  seme 
optimum  frequency  for  transmission  of  sound  in  a  channel  of 
given  thickness,  lying  somewhere  above  the  cutoff  frequency 
and  below  the  frequency  range  where  absorption  and  scat¬ 
tering  losses  become  large. 


absorption.  The  degree  of  attenuation  is  strongly  dependent 
on  the  frequency  of  the  sound,  with  higher  frequencies 
experiencing  greater  attenuation.  Thorp  [Bef.  4]  developed 
extensive  tables  to  describe  the  attenuation  of  sound  as  a 
function  of  frequency. 

Sound  can  be  scattered  from  rough  boundaries  at  the  sea 
surface  or  bottom,  or  from  discontinuities  at  the  edges  of 
the  layers.  Bubbles,  marine  organisms  and  suspended  partic¬ 
ulate  matter  will  also  scatter  sound  [Bef.  3].  Scattering 
is  generally  more  important  for  higher  frequencies  than  for 
lower  ones. 

Sound  energy  can  also  be  lost  through  absorption. 
Viscosity  effects  in  the  interaction  of  water  molecules  as 
the  sound  wave  passes  account  for  one  absorption  mechanism. 
Chemical  reactions  are  also  important  sources  of  absorptive 
loss.  Of  the  many  chemical  constituents  of  seawater, 
certain  compounds  appear  to  undergo  continuous  dissociation 
into  ions  and  reassociation.  This  relaxation  mechanism 
takes  up  acoustic  energy.  The  most  important  absorption 
mechanism  for  frequencies  near  1000  Hz  appears  to  be  the 
borate-boric  acid  relaxation  mechanism  [Bef.  5]. 

In  general,  attenuation  increases  logarithmically  with 
frequency.  The  attenuation  losses  due  to  such  mechanisms 
are  so  small  at  low  frequencies  that  the  measurements  must 
be  made  over  extremely  long  ranges  to  detect  any  losses. 
The  precision  with  which  the  measurements  must  be  made  and 
the  inhomogeneities  in  the  ocean  further  complicate  the 
task.  Orick  [Bef.  2]  provides  an  excellent  discussion  of 
absorption  and  attenuation. 

Two  factors  are  operating  with  opposing  effects  on  sound 
being  propagated  through  sound  channels.  First,  attenuation 
mechanisms  tend  to  cause  losses  of  acoustic  energy.  Second, 
the  focusing  effect  of  the  channel  tends  to  reduce  the 
losses.  In  the  ideal  case,  a  sound  channel  bounded  above 
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formula  is  given  by  Kinsler,  Frey, 

[Ref.  3] 

p  ~  2x  I05  /  3 

hCO  ~  /  D  '2 


Coppens  and  Sanders 


(egn  2.2) 


for  0  in  aeters.  This  yields  a  cutoff  frequency  of  about 
280  Hz  for  this  channel.  All  of  these  calculations  agree 
within  useful  ranges,  since  the  cutoff  frequency  does  not 
mark  an  abrupt  end  to  propagation. 


C.  TBAHSHISSIOH  LOSS 

Several  oechanisas  act  to  decrease  the  energy  of  an 
acoustic  wave  as  the  wavefront  moves  farther  away  from  its 
source.  The  most  obvious  loss  mechanism  is  spreading.  Once 
an  outwardly  propagating  wavefront  reaches  a  reasonable 
distance  from  a  sufficiently  small  source,  the  wavefront  can 
be  considered  for  all  practical  purposes  to  be  spherical. 
Thus,  the  energy  in  the  wave  is  spread  out  over  an  area  that 
increases  as  the  square  of  the  range  from  the  source.  In 
decibels,  spherical  spreading  alone  would  produce  a  trans¬ 
mission  loss  (TL)  of 


TLspher  Io9ioR 


(eqn  2.3) 


If  the  sound  is  somehow  contained  so  that  spreading  may 
only  take  place  horizontally,  as  is  the  case  when  the  sound 
source  is  inside  a  channel  or  duct,  the  spreading  is  cylin¬ 
drical  rather  than  spherical  and  the  loss  is  given  by 


TLcyl  -  10  iog|0 

Other  mechanisms  cause  additional  energy  to  be  lost. 
These  attenuation  mechanisms  include  scattering  and 


(eqn  2.4) 
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B.  ACOUSTIC  WAVEGUIDES 


Sound  channels  act  as  acoustic  waveguides.  The  ray 
theory  of  sod'nd  propagation  predicts  that  sound  can  be 
refracted  and  trapped  in  ducts.  Trapping  depends  on  the 
thickness  of  the  duct  and  the  wavelength  of  the  sound;  if 
the  wavelength  is  "too  large"  to  be  contained  within  the 
vertical  dimension  of  the  duct,  the  sound  will  propagate 
poorly.  Frequencies  corresponding  to  those  longer  wavel¬ 
engths  will  propagate  poorly,  if  at  all;  higher  frequencies 
with  shorter  wavelengths  will  propagate  much  better.  Normal 
■ode  theory  provides  a  tool  for  understanding  these  propaga¬ 
tion  effects  £Bef.  3].  The  transition  from  poor  to  good 
propagation  is  marked  by  the  cutoff  frequency.  This 
frequency  does  not  mark  an  abrupt  end  to  propagation  below 
which  sound  cannot  be  transmitted  at  all;  rather  it  is  a 
lower  boundary  below  which  propagation  is  less  effective. 

Acoustic  propagation  in  a  sound  channel  is  exactly  anal¬ 
ogous  to  the  propagation  of  radio  waves  in  a  ground-based 
duct.  Urick  [Eef.  2]  provides  an  equation  from  the  theory 
of  radio  propagation  for  the  maximum  wavelength  which  will 
"fit"  in  such  a  duct.  Dosso  and  Chapman  [Bef.  1]  use  this 
equation  to  calculate  the  cutoff  frequency  for  their 
80-meter  thick  shallow  sound  channel  as  about  250  Hz. 

Urick  gives  another  equation  derived  from  radio  propaga¬ 
tion  theory  to  calculate  the  cutoff  frequency;  [Bef.  2] 

Amox  =4.7x|03  H  /2  (eqn  2.1) 

where  ^max  is  the  oaximum  wavelength  of  the  trapped  sound 
in  feet,  and  H  is  the  thickness  of  the  duct  in  feet. 
Converted  into  metric  units  and  using  a  sound  velocity  of 
1483  m/s,  this  formula  yields  a  cutoff  frequency  of  about 
243  Hz  for  Dosso  and  Chapman’s  80-meter  channel.  A  third 


A  typical  sound  velocity  profile  like  that  of  Figure  2.1 
[Bef.  3]  shows  a  surface  layer  of  seawater  which  is  subject 
to  the  actions,  of  sun,  wind  and  precipitation.  This  layer 
is  often  well-aired,  but  its  temperature,  salinity  and 
thickness  may  vary  greatly  from  place  to  place  and  from  day 
to  day.  Typically  the  surface  layer  shows  a  constant  or 
slightly  increasing  sound  velocity  profile.  The  surface 
layer  may  extend  to  a  depth  of  several  tens  of  meters.  At 
some  point  below  the  surface  layer  a  strong  thermal  gradient 
may  show  up.  Temperature  and  therefore  sound  velocity 
decrease  sharply  with  depth,  creating  a  seasonal  thermo- 
cline.  Deeper  still  lies  the  main  thermocline,  formed  by 
another  strong  negative  thermal  gradient.  Below  this  is  the 
deep  isothermal  layer,  which  extends  to  the  sea  floor.  In 
this  layer  salinity  is  essentially  constant,  as  is  tempera¬ 
ture.  Pressure  and  velocity  increase  very  nearly  linearly 
with  depth.  The  layers  bounded  above  by  the  main  thermo¬ 
cline  and  below  by  the  increasing  velocity  gradient  are 
called  the  deep  sound  channel  {DSC)  .  Sound  can  be  trapped 
in  this  channel  and  refracted  repeatedly  toward  the  velocity 
minimum,  which  is  called  the  DSC  axis.  This  axis  may  occur 
at  depths  ranging  from  several  hundred  to  over  1000  meters. 

Certain  environmental  conditions  can  cause  a  sound 
velocity  profile  like  that  shown  in  Figure  2.2.  A  combina¬ 
tion  of  winds,  high  precipitation  and  low  evaporation  such 
as  is  found  in  the  northeastern  Pacific  Ocean  can  lead  to 
the  formation  of  a  secondary  or  shallow  sound  channel  (SSC) 
just  below  the  seasonal  thermocline.  Just  as  sound  can  be 
trapped  in  the  DSC,  sound  can  also  be  trapped  in  the  SSC 
[Ref.  1]. 


A.  SCU1D  II  THE  SEA 


The  Banner  in  which  sound  is  transaitted  through  a 
aediua  depends  on  the  velocity  with  which  it  propagates. 
Because  of  the  greater  density  of  seawater,  sound  is  trans¬ 
aitted  auch  faster  in  the  ocean  than  in  the  air.  The 
velocity  of  sound  at  any  given  point  in  the  ocean  is  deter- 
ained  by  the  temperature,  salinity  and  pressure.  Since 
pressure  varies  linearly  with  water  depth,  these  two  quanti¬ 
ties  can  be  treated  as  interchangeable.  Teaperature  and 
depth  are  generally  the  doainant  factors  affecting  the  speed 
of  sound.  Thus  the  velocity  of  sound  in  the  ocean  varies  in 
all  diaensions:  horizontally,  vertically  and  teaporally 
[Ref.  2]. 

The  sound  velocity  can  be  aeasured  directly  using  a 
velociaeter  or  indirectly  using  a  bathythermograph  (BT)  to 
aeasure  temperature  as  a  function  of  depth.  Numerous 
researchers  have  provided  eapirical  formulas  for  calculating 
velocity  froa  teaperature,  salinity  and  depth  [Ref.  3].  In 
either  case  the  result  is  a  record  of  the  sound  velocity  as 
it  changes  with  depth,  or  sound  velocity  profile  (SVP)  . 

Snell’s  Law  of  acoustics  is  exactly  analogous  to  Snell’s 
Law  of  optics.  When  a  ray  of  light  or  sound  passes  froa  one 
aediua  into  another,  the  ray  is  bent  or  refracted  toward  the 
aediua  in  which  it  travels  at  the  slower  speed.  The  ocean 
may  be  thought  of  as  composed  of  aany  layers  of  water,  in 
each  of  which  sound  has  a  different  velocity.  Where  the  STP 
shows  a  velocity  minimum,  sound  is  refracted  toward  that 
depth,  resulting  in  a  focusing  of  energy  at  that  velocity 
ainiaua.  This  depth  is  often  referred  to  as  a  sound  channel 
axis . 


Dosso  and  Chapman  were  able  to  obtain  good  agreement  between 
their  experimental  results  and  the  predictions  of  the  model. 
A  full  discussion  of  Dosso  and  Chapman's  experiment  is 
contained  in  Chapter  III. 

This  study  examines  the  Dosso  and  Chapman  experiment, 
then  attempts  to  reproduce  their  results  using  the  Naval 
Postgraduate  School's  PE  and  FACT  9H  models.  Chapter  IV 
compares  these  predictions  with  the  Dosso  and  Chapman 
results,  supporting  the  use  of  the  two  models  to  predict  the 
optimum  frequency  for  sound  propagation  within  a  shallow 
sound  channel.  The  models  are  then  applied  to  two  sound 
velocity  profiles  of  interest  to  the  Naval  Postgraduate 
School's  Environmental  Acoustic  Research  Group.  These 
profiles  were  obtained  under  summer  conditions  in  a  deep 
ocean  basin  (4000  meters)  at  high  latitude.  Chapter  V 
describes  the  results  of  the  NPS  PE  model  as  applied  to 
these  profiles.  Chapter  VI  discusses  the  FACT  9H  model  and 
its  predictions  for  the  Dosso  and  Chapman  profiles  and  the 
high- latitude  summer  profiles.  Chapter  VII  then  summarizes 
this  study's  results  and  conclusions  and  presents  recommen¬ 
dations  for  further  investigation. 


The  first  model,  the  split-step  Parabolic  Equation  (PE) 
model,  uses  a  parabolic  approximation  to  the  basic  wave 
equation  and  solves  for  pressure  levels  as  a  function  of 
range.  It  is  considered  to  be  a  low-frequency  (below  800 
Hz)  model,  since  higher  frequencies  require  large  amounts  of 
computational  time.  The  PE  model,  a  full-wave  solution,  is 
considered  a  benchmark  against  which  other  transmission  loss 
models  are  tested.  It  is  not  commonly  used  in  support  of 
fleet  operations  due  to  the  large  computational  times 
required. 

The  second  model,  the  Fast  Asymptotic  Coherent 
Transmission  loss  model  version  9H  (FACT  9H),  employs  ray¬ 
tracing  techniques  to  predict  the  amount  of  acoustic  energy 
arriving  at  a  receiver  as  a  function  of  range.  This  model 
is  less  rigorous  in  predicting  the  effects  of  different 
types  of  boundaries  and  it  does  not  allow  for  variations  of 
the  sound  velocity  profile  with  range.  It  is  simpler  and 
requires  less  computational  time.  FACT  is  generally  consid¬ 
ered  more  appropriate  than  the  PE  model  for  frequencies 
above  800  Hz.  The  FACT  model  is  widely  used  in  acoustic 
transmission  loss  predictions  provided  to  fleet  users  in  the 
0.  S.  Navy.  Both  the  PE  model  and  FACT  9H  are  available  at 
the  Naval  Postgraduate  School,  accessed  through  the  IBM  3033 
computer. 

In  1984  Dosso  and  Chapman  [Bef.  1]  reported  an  experi¬ 
ment  on  propagation  loss  that  they  had  conducted  in  the 
northeast  Pacific  Ocean.  Sound  velocity  profiles  indicated 
the  presence  of  a  shallow  sound  channel  overlying  a  deep 
sound  channel.  They  placed  sources  of  known  magnitude  in 
the  shallow  sound  channel,  with  receivers  in  the  shallow 
sound  channel  and  the  deep  sound  channel.  They  measured 
propagation  losses  over  ranges  extending  from  zero  to  35  km 
for  propagation  both  within  the  SSC  and  across  the  boundary 
with  the  DSC.  Using  a  split-step  parabolic  equation  model. 


resulting  from  winds  and  turbulence  which  may  extend  from 
the  surface  to  depths  of  several  tens  of  meters.  Below  the 
mixed  layer  a-  shallow  sound  channel  (SSC)  is  formed  by  a 
combination  of  temperature  and  salinity  gradients,  resulting 
in  a  sound  speed  minimum  (channel  axis)  at  roughly  80  to  200 
meters.  The  lower  edge  of  the  SSC  is  bounded  by  an  area  of 
rapidly  increasing  sound  speed.  Below  the  SSC  the  combined 
effects  of  decreasing  temperature  and  increasing  pressure 
can  also  result  in  a  second  or  deep  sound  channel  (DSC)  with 
an  axis  at  about  800  to  1000  meters  in  mid-latitudes.  Sound 
channels  act  as  acoustic  waveguides,  trapping  certain 
frequencies  and  allowing  them  to  propagate  especially  well, 
while  other  frequencies  "leak”  out  of  the  channel  and  propa¬ 
gate  poorly.  A  more  complete  discussion  of  the  theoretical 
background  for  this  paper  is  contained  in  Chapter  II. 

The  various  layers  and  channels  discussed  above  occupy 
the  top  several  hundred  to  one  thousand  meters  of  the 
oceans.  Submarines  operate  within  these  depths.  These 
features  of  the  sound  velocity  profile  strongly  affect  the 
propagation  of  the  submarine's  acoustic  signature,  and  the 
active  and  passive  operation  of  the  submarine's  sonar  equip¬ 
ment  and  that  of  its  hunters.  Antisubmarine  warfare  (ASW) 
operators  must  understand  these  effects  and  have  some  means 
of  predicting  them.  With  reliable  predictions  operators 
will  be  able  to  select  frequencies  and  modes  of  operation  to 
maximize  the  effectiveness  of  their  sonar  equipment  in 
detecting  and  tracking  the  submarine. 

B.  OBJECTIVE 

This  study  seeks  to  use  two  acoustic  transmission  loss 
models  to  predict  the  behavior  of  sound  in  typical  shallow 
sound  channels.  Both  models  are  available  at  the  Naval 
Postgraduate  School  (NPS) ,  accessed  through  the  IBM  3033 
computer. 
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I.  IBTBODOCTIOH 


1.  BACKGROUND 

Despite  technological  advances  in  ship  silencing,  the 
acoustic  signature  produced  by  a  submarine’s  main  propulsion 
system  and  auxiliary  machinery  remains  the  single  most 
easily  exploited  tool  for  detecting,  tracking  and  localizing 
the  submarine.  The  transmission  of  sound  through  the  ocean 
is  predictable  given  a  sufficiently  well-defined  knowledge 
of  the  environment,  the  source  and  the  receiver.  The  trans¬ 
mission  of  sound  in  the  sea  depends  on  the  sound  velocity 
profile  (SVP)  .  Sound  velocity  in  turn  depends  on  the 
temperature,  salinity  and  depth  of  the  water.  Transmission 
phenomena  can  become  extremely  complex,  even  for  relatively 
simple  sound  velocity  profiles. 

Acoustic  transmission  loss  models  have  been  developed 
for  implementation  on  digital  computers.  The  models  gener¬ 
ally  take  one  of  two  approaches:  either  solving  the 
acoustic  wave  eguaticn  or  employing  ray-tracing  techniques. 
Both  approaches  require  assumptions,  simplifications  and 
approximations.  For  real  situations  some  models  will 
predict  transmission  loss  more  accurately  than  others 
depending  on  the  computational  technique  involved,  the 
method  of  solution  and  the  extent  to  which  the  model  accepts 
real  conditions  for  boundaries  at  the  sea  surface,  the  sea 
floor  and  between  different  water  masses. 

The  situation  of  interest  here  is  one  of  the  more 
complicated  cases:  the  shallow  sound  channel.  This  situ¬ 
ation  develops  in  higher  latitudes  of  the  world's  oceans, 
where  precipitation  is  heavy  and  evaporation  is  small.  A 
sound  velocity  profile  reveals  the  presence  of  a  mixed  layer 
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energy  levels.  Bottom  interactions  could  not  be  isolated 
from  aore  distant  signals.  Propagation  loss  was  calculated 
as  known  source  level  minus  the  received  level  measured  at 
the  hydrophones. 

B.  RESULTS  AID  DISCDSSION 

For  the  shallow  sound  channel  [Ref.  1  ],  measured  propa¬ 
gation  loss  was  plotted  as  a  function  of  range  for  the  1/3 
octave  hands.  The  propagation  loss  demonstrated  strong 
dependence  on  range,  as  expected.  An  optimum  frequency  for 
propagation  of  about  800  Hz  was  observed. 

Figure  3.4  shows  the  propagation  loss  measured  by  Dosso 
and  Chapman  for  both  source  and  receiver  located  within  the 
shallow  sound  channel.  The  solid  curve  represents  geome¬ 
trical  spreading  loss,  calculated  from 

H  -  1 0  lo0lo  ( R  •  R0 )  (eqn  3*1) 

where  range  R  and  ocean  depth  R0  are  measured  in  meters. 
Figure  3.4  (a)  represents  propagation  loss  below  the  optimum 
frequency  cf  800  Hz;  Figure  3.4  (b)  represents  propagation 
loss  above  800  Hz.  The  close  agreement  between  the  geome¬ 
trical  loss  curve  and  the  800  Hz  data  over  much  of  the  range 
is  noteworthy. 

Figure  3.5  displays  the  propagation  loss  within  the 
shallow  sound  channel  for  the  optimum  propagation  freguency 
of  800  Hz,  with  the  effects  of  chemical  absorption  removed. 
Again  the  solid  line  represents  geometrical  spreading  for 
comparison.  Note  the  abrupt  increase  in  losses  over  the 
8-20  km  range.  Over  this  range  the  sound  channel  axis  deep¬ 
ened,  so  that  the  sources  were  no  longer  within  the  duct. 

Figure  3.6  presents  the  same  data  plotted  as  a  function 
of  frequency  for  three-point  range  averages  centered  at  24 


and  33  km.  The  solid  circles  represent  actual  measurements; 
the  open  circles  represent  values  adjusted  upward  to  remove 
the  effects  of  chemical  absorption.  The  dotted  lines  repre¬ 
sent  the  expected  less  from  geometrical  spreading  alone  at 
800  Hz.  This  display  shows  more  clearly  the  optimum 
frequency  of  800  Hz  and  how  closely  the  losses  for  this 
frequency  approximate  the  losses  due  to  geometrical 
spreading  alone.  It  was  clear  that  frequencies  above  and 
below  800  Hz  experienced  greater  losses.  Diffraction  or 
"leakage”  from  the  shallow  sound  channel  accounted  for  the 
losses  at  lower  frequencies.  High  frequencies  experienced 
losses  due  to  chemical  absorption  mechanisms  and  scattering. 

C.  COMPARISON  BETWEEN  SSC  AMD  DSC  PROPAGATION 

The  deep  sound  channel  trapped  lower  frequencies  (longer 
wavelengths)  than  the  shallow  sound  channel  because  the  DSC 
is  much  thicker.  Figure  3.7  demonstrates  the  effects  of 
receiver  depth  on  propagation  loss.  The  sources  were 
located  at  98  meters,  near  the  shallow  sound  channel  axis. 
One  receiver  was  also  near  the  shallow  sound  channel  axis  at 
101  meters.  The  other  receiver  was  positioned  at  417 
meters,  near  the  deep  sound  channel  axis.  At  160  Hz  sound 
leaked  easily  out  of  the  shallow  sound  channel  but  was  well 
trapped  within  the  deep  sound  channel.  At  315  Hz,  just 
above  the  cutoff  freguencj  for  the  shallow  sound  channel, 
propagation  in  the  channel  improved.  At  630  Hz  sound  was 
effectively  trapped  in  the  shallow  sound  channel,  and  propa¬ 
gation  loss  at  the  deep  receiver  was  significantly  greater 
than  at  the  shallow  receiver. 

Dosso  and  Chapman  plotted  relative  gain  between  the  two 
receivers  according  to  the  relationship 


Relative  Gain  nTI  — T| 

^)EEP  1  '-SHALLOW 


(eqn  3.2) 


Figure  3.8  displays  the  results  for  a  four-point  range 
average  centered  at  24  km.  The  cutoff  frequency  of  250  Hz 
is  clearly  demonstrated.  Below  250  Hz  the  SSC  fails  to  trap 
the  sound  effectively,  giving  larger  transmission  losses 
than  the  DSC  and  a  negative  relative  gain.  Above  250  Hz  the 
SSC  traps  the  sound  tetter,  producing  smaller  transmission 
losses  than  the  DSC  and  a  positive  relative  gain.  Larger 
propagation  losses  in  the  SSC  below  25  Hz  result  from 
surface  decoupling  effects. 

D.  HODELIBG  PEOPAGATIGH  LOSS  1H  THE  SSC 

Geometrical  spreading  is  a  simple  model,  but  it  is  not 
able  to  predict  the  frequency  and  range  dependence  that 
Dosso  and  Chapman  observed.  A  split-step  parabolic  equation 
(PE)  model  was  used  to  calculate  propagation  loss  based  on 
the  sound  velocity  profiles  measured  at  ranges  of  zero,  16 
and  33  km  along  the  track.  [Ref.  1].  Attempts  to  model 
propagation  using  only  these  three  profiles,  however,  were 
unsatisfactory.  The  16  km  profile  indicated  that  the  axis 
of  the  SSC  had  deepened,  so  that  the  sources  were  no  longer 
contained  in  the  channel.  This  is  reflected  in  the  increased 
propagation  loss  over  the  5-20  km  range  in  Figure  3.5.  The 
improved  propagation  at  and  beyond  20  km  indicated  to  Dosso 
and  Chapman  that  the  axis  of  the  SSC  must  have  risen  to 
about  100  meters  again.  For  this  reason  they  used  the 
profile  recorded  at  33  km  for  the  20  km  range  also.  Dosso 
and  Chapman  also  tzuncated  their  four  sound  velocity 
profiles  at  a  depth  of  600  meters,  introducing  an  acousti¬ 
cally  transparent  bottom  at  that  depth  [Ref.  7].  This  step 
had  the  practical  effects  of  limiting  the  Fourier  transform 
matrix  to  a  manageable  size  and  shortening  computational 
time  somewhat.  It  may  also  have  had  some  effect  on  the 
predicted  losses. 
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With  these  four  sound  velocity  profiles  Dosso  and 
Chapman's  aodel  predicted  losses  vith  an  acceptable  degree 
of  accuracy  at  25  Hz  and  630  Hz  (the  practical  lower  and 
upper  linits  of  their  model) .  But  near  the  cutoff  frequency 
of  250  Hz,  the  predicted  losses  were  not  in  such  good  agree¬ 
ment  with  observed  losses.  The  model  did,  however,  predict 
strong  dependence  on  both  frequency  and  range  as  Ocssc  and 
Chapman  observed.  Figure  3.9  shows  the  model  results 
compared  to  actual  measurements. 

B.  SOHBABT  OF  DOSSO  AHD  CHAPHAH'S  EXPEBIMEHT 

Dosso  and  Chapman  point  out  that  the  shallow  sound 
channel  behaved  like  an  acoustic  waveguide  with  optimum 
propagation  at  800  Hz,  with  increased  losses  at  lower 
frequencies  due  to  diffraction  and  at  higher  frequencies  due 
to  absorption  and  scattering  £Bef.  1].  The  propagation  was 
found  to  be  very  sensitive  to  changes  in  the  environment. 
The  PE  model  that  Dosso  and  Chapman  used  was  able  to  account 
for  these  range-dependent  effects  with  good  agreement  with 
experimental  results  at  low  and  near-optimum  frequencies. 
At  intermediate  frequencies  the  model  showed  range  and 
frequency  dependent  effects,  but  the  predicted  losses  were 
greater  than  experimental  results  indicated.  Dosso  and 
Chapman  emphasize  that  more  accurate  modeling  of  propagation 
losses  requires  a  much  finer  sampling  of  the  environment; 
that  is,  sound  velocity  profiles  at  much  closer  intervals 
than  16  kilometers. 
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Figure  3.1  Dosso  and  Chapaan’s  Sound  Telocity  Profile 
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Figure  3.7  Propagation  Loss  vs  Range  for  Deep 
and  Shallow  Receivers  (Dosso  and  Chapman) • 
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Fiaure  3.8  Relative  Gain  between  Deep  and  Shallow  Receivers 

(Dosso  and  Chapaanj . 
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I L.  THE  HES  P ABA BO IIC  EQUATION  MODEL 


Dosso  and  Chapman  reported  good  results  [Ref.  1]  using 
their  split-step  Parabolic  Equation  (PE)  model  to  predict 
propagation  loss  in  a  range-dependent  environment  where  both 
a  shallow  sound  channel  and  a  deep  sound  channel  existed. 
The  first  step  in  this  study  was  to  attempt  to  reproduce 
those  results  using  the  PE  model  available  at  Naval 
Postgraduate  School  Mcnterey. 

A.  MODEL  DESCRIPTION 

The  split-step  Parabolic  Equation  (PE)  acoustic  trans¬ 
mission  loss  model  was  chosen  for  study  because  it  produces 
a  numerical  solution  for  an  approximation  to  the  wave  equa¬ 
tion  and  is  generally  considered  to  provide  the  most  accu¬ 
rate  tool  for  predicting  transmission  loss.  It  is  often 
used  as  a  benchmark  against  which  to  compare  the  performance 
of  other  transmission  loss  models.  It  can  deal  with  range- 
dependent  environments  (multiple  sound  velocity  profiles) 
and  with  some  interaction  of  sound  with  the  sea  floor.  Its 
drawbacks  are  its  relatively  low-frequency  capability, 
inability  to  handle  all  bottom  loss  conditions,  and  lack  of 
flexibility  in  dealing  with  boundary  conditions  at  the  sea 
surface  (the  model  uses  a  flat  pressure-release  boundary 
which  does  not  account  for  the  effects  of  wind  and  waves) . 
This  PE  model  is  implemented  on  the  IBM  3033  computer  at  the 
Naval  Postgraduate  School  (NPS)  and  is  accessed  through  an 
interactive  program  developed  and  managed  by  the  school’s 
oceanography  department. 

The  FORTRAN  source  code  for  this  PE  model  was  written  by 
the  Acoustic  Environmental  Support  Detachment  (A ESD)  of  the 
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Figure  4.5  Transmission  Loss  vs  Bange,  25  Hz  (BPS  PE  Model) . 
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Figure  4.4  Transuission  Loss  rs  Range ,  100  Hz  (NPS  PE  Model) . 
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Figure  9.1  Transmission  loss  vs  Range,  2000  Hz 

(NPS  PE  Hodei) . 
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demonstrated  an  appropriate  cutoff  frequency  around  250  Hz, 
tut  with  unexpected  peaks  at  frequencies  well  below  cutoff. 

Seme  of  these  anomalous  results  from  the  NPS  PE  model 
■ay  arise  froa  using  full-depth  sound  velocity  profiles. 
Seasoning  froa  ray-tracing  theory,  a  fully  absorbing  bottom 
placed  at  2500-2900  meters  aay  have  allowed  sufficient  depth 
for  soae  deep  propagation  paths  to  be  refracted  back  toward 
the  receiver.  The  fully  absorbing  bottom  at  600  aeters 
which  Dosso  and  Chapaan  used  in  their  aodeling  aay  have 
eliminated  propagation  by  these  paths.  Interference  effects 
between  these  paths  and  the  ducted  propagation  may  account 
for  seme  of  the  differences  between  Dosso  and  Chapman's 
results  and  the  NPS  PE  results. 
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about  18  ka,  followed  by  better  propagation  at  longer 
ranges.  Their  predictions  correspond  rather  well  to  the 
aeasureaents. _  On  the  other  hand.  Figure  4.12  predicts 
several  dB  better  propagation  for  630  Hz  than  was  observed 
for  ranges  less  than  15  ka  and  several  dB  worse  than 
observed  at  ranges  beyond  20  ka.  Again  the  compensation 
using  the  20  km  profile  was  less  successful  for  the  NPS  PE 
aodel  than  for  Oosso  and  Chapman's  model.  The  steep  drop 
and  failure  to  recover  in  the  transmission  loss  curves  are 
even  more  apparent  in  Figure  4.13  for  400  Hz  and  Figure  4.14 
for  250  Hz.  In  Figure  4.15  for  25  Hz  the  observed  values 
Batch  quite  well  with  predicted  values  out  to  about  15  km. 
At  that  point  the  prediction  goes  wrong  again  and  predicted 
values  fall  10-15  dB  below  measured  values. 

C.  DISCUSSION  OF  BPS  PE  RESULTS 

Discounting  the  multipath  effects  on  transmission  loss 
which  produce  extreme  fluctuations  in  transmission  loss  with 
both  range  and  frequency,  it  appears  valid  to  use  the  NPS 
Parabolic  Equation  model  to  estimate  optimum  frequency  for 
sound  propagation  in  a  shallow  sound  channel.  The  plot  of 
transmission  loss  versus  frequency  (Figure  4.7)  shows  an 
optimum  frequency  range  of  500-700  Hz.  This  is  acceptable 
in  comparison  with  Cosso  and  Chapman's  measured  optimum 
frequency  of  800  Hz  for  the  same  set  of  sound  velocity 
profiles.  This  is  also  good  enough  for  active  and  passive 
sonar  operators  to  use  in  selecting  the  best  frequency 
ranges  for  their  equipment. 

The  NPS  PE  aodel  was  not,  however,  a  good  estimator  for 
absolute  levels  of  transmission  loss.  The  NPS  PE  curves 
showed  large  differences  in  transmission  loss  levels  when 
compared  to  Dosso  and  Chapman's  measurements.  For  relative 
gain  between  deep  and  shallow  receivers  the  NPS  PE  model 


trapping  and  leaking  are  better  balanced.  There  is  still  a 
10  dB  better  propagation  for  the  15-17.5  km  range  when  both 
source  and  receiver  are  in  the  channel.  Figure  4.10  demon¬ 
strates  significant  transmission  loss  for  160  Hz,  which  is 
well  below  cutoff.  In  this  case  the  deep  receiver  experi¬ 
ences  better  propagation  by  5-10  dB  over  almost  the  entire 
range,  indicating  that  the  shallow  sound  channel  failed  to 
trap  this  freguency  effectively. 

It  should  especially  be  noted  in  Figures  4.8  and  4.9 
that  ancmalous  propagation  is  occurring  in  the  17-25  km 
range.  This  may  be  due  to  the  change  in  depth  of  the 
shallow  sound  channel  that  Dosso  and  Chapman  noted  and  tried 
to  correct  by  adding  the  sound  velocity  profile  at  20  km. 
This  compensation  was  apparently  much  less  successful  with 
the  NFS  EE  model  than  with  Dosso  and  Chapman's  model. 

Figure  4.11  corresponds  to  Figure  3.8.  Relative  gain  is 
calculated  according  to  equation  3.2,  following  Dosso  and 
Chapman's  procedure.  Their  results  were  good  for  24  km, 
showing  a  sudden  improvement  above  the  250  Hz  cutoff 
frequency  when  both  source  and  receiver  were  in  the  SSC. 
Ihe  curve  for  the  NPS  PE  model  at  the  24  km  range  does  not 
produce  good  results,  probably  due  to  the  failure  of  the 
added  profile  at  20  km  to  compensate  for  the  changing  depth 
of  the  SSC.  The  15  km  curve,  however,  is  taken  before  the 
SSC  deepened  and  indicates  better  results.  The  250  Hz 
cutoff  frequency  is  clearly  indicated.  The  NPS  results, 
however,  show  unexpected  peaks  at  40-50  Hz  and  at  100  Hz, 
where  the  shallow  sound  channel  receiver  should  not  be  doing 
as  well  as  the  deep  receiver.  Other  frequencies  below 
cutoff  behave  as  expected. 

Together  Figures  4.12  through  4.15  correspond  to  Figure 
3.9,  where  Dosso  and  Chapman  plotted  their  PE  model 
predicticns  along  with  experimentally  measured  data.  Note 
that  Dosso  and  Chapman's  630  Hz  curve  shows  a  minimum  at 


The  NFS  model  allows  the  user  to  select  whether  attenuation 
will  be  considered.  Figure  4.6  represents  the  NPS  PE  iodel 
prediction  for.  800  Hz  without  attenuation.  Again  the  steep 
drop  in  transmission  occurs  between  15  and  20  km.  A  compar¬ 
ison  with  Figure  4. 2  shows  that  the  two  curves  are  essen¬ 
tially  similar,  but  attenuation  accounts  for  a  1-2  dB 
greater  loss  in  Figure  4.2. 

Ocsso  and  Chapman* s  experimental  data  yielded  good 
results  when  propagation  loss  was  plotted  as  a  function  of 
frequency  for  24  and  33  km  (Figure  3.6).  Recall  that  Dosso 
and  Chapman  measured  1/3  octave  band  levels  and  then  calcu¬ 
lated  three-point  range  averages  to  obtain  each  curve.  The 
NPS  PE  model  produces  transmission  loss  predictions  for 
discrete  frequencies.  To  obtain  transmission  loss  for  the 
desired  ranges,  values  were  read  from  the  computer-generated 
PE  model  graphic  output  at  the  appropriate  ranges  and  at  one 
kilometer  on  either  side.  The  three  values  for  each  range 
were  averaged  to  produce  the  values  used  in  plotting  Figure 
4.7.  Again  multipath  interference  effects  can  account  for 
the  fluctuations  with  range  and  frequency.  While  an  optimum 
frequency  of  800  Hz  is  easy  to  read  from  Figure  3.6,  it  is 
more  difficult  to  state  an  optimum  frequency  from  Figure 
4.7.  Frequencies  in  the  500-700  Hz  range  seem  to  have 
propagated  best. 

Dosso  and  Chapman  compared  results  from  shallow  sound 
channel  and  deep  sound  channel  receivers  in  Figure  3.7. 
Figures  4.8  through  4.10  present  the  corresponding  results 
from  the  NPS  PE  model.  As  expected  for  a  frequency  well 
above  cutoff,  the  630  Hz  curves  (Figure  4.8)  show  better 
propagation  where  both  source  and  receiver  were  within  the 
shallow  sound  channel,  but  much  greater  loss  (up  to  10-20 
dB)  where  the  receiver  was  in  the  deep  sound  channel.  The 
two  curves  are  in  closer  agreement  for  315  Hz  (Figure  4.9) 
as  would  be  expected  for  a  frequency  near  cutoff,  where 


Ccnpare  Figures  4.1  through  4.5  with  Figure  3.4.  The 
smooth,  downward  trending  curve  represents  transmission  loss 
due  to  geometric  spreading  alone,  calculated  from  eguation 
3.1.  The  measured  1/3  octave  band  levels  of  Figure  3.4  show 
signs  of  smoothing  because  many  frequencies  have  been  added 
together,  lessening  the  multipath  interference  effects  that 
become  apparent  for  discrete  frequencies  [Bef.  1]. 

It  is  apparent,  however,  that  frequencies  below  the 
cutoff  frequency  of  250  Hz  propagate  poorly  within  the 
shallow  sound  channel  and  that  higher  frequencies  propagate 
better.  in  interesting  outcome  here  is  that  the  2000  Hz 
experiences  less  transmission  loss  than  the  experimentally 
determined  optimum  frequency  of  800  Hz  at  ranges  greater 
than  17.5  km.  Note  also  that  the  three  highest  frequencies 
all  show  sudden  steep  increases  in  transmission  loss  in  the 
17-20  km  range  and  that  the  levels  never  really  recover. 
This  is  the  range  at  which  Oosso  and  Chapman  noted  a  change 
in  their  sound  speed  profiles  so  that  their  sources  were  no 
longer  within  the  shallow  sound  channel. 

Dosso  and  Chapman  observed  an  increase  in  transmission 
loss  over  roughly  the  same  ranges  (Figure  3.4),  folloved  by 
better  propagation  at  ranges  beyond  20  km.  They  were  able 
to  model  these  results  successfully  by  adding  another  sound 
velocity  profile  at  20  km,  using  the  same  profile  they  had 
measured  at  33  km.  lith  the  same  input  parameters  the  NPS 
PE  model  yielded  a  much  steeper  drop  and  did  not  indicate  a 
return  to  levels  expected  from  the  experimental  data. 

Oosso  and  Chapman’s  Figure  3.5  compares  to  Figure  4.6. 
For  Figure  3.5  Oosso  and  Chapman  plotted  measured  data,  with 
losses  due  to  chemical  absorption  removed.  This  figure 
clearly  shows  the  close  correspondence  between  the 
measurements  for  the  optimum  frequency  of  800  Hz  and  the 
losses  predicted  for  geometrical  spreading  alone.  The  only 
anomalies  arise  over  the  8-20  ’  **  range,  as  discussed  above. 


solved  numerically  using  the  Tappert-Hardin  split-step 
Fourier  algorithm  tc  yield  a  solution  for  the  entire 
acoustic  pressure  field  as  a  function  of  range  and  depth 
[Ref.  6]. 

The  ?E  solution  is  best  where  the  source  and  receiver 
are  arranged  so  that  the  path  between  the  two  makes  a  small 
angle  with  the  hcrizcntal.  At  larger  angles,  the  group  and 
phase  velocities  of  the  computed  solution  become  subject  to 
errors  and  produce  inaccurate  representations  in  the  inter¬ 
ference  pattern.  This  affects  the  accuracy  of  the  predicted 
transmission  loss  for  a  given  range.  Computing  time 
increases  quickly  with  frequency.  The  numerical  solution  is 
limited  by  the  practical  Fourier  transform  size,  which  means 
that  the  input  environment  governs  the  extent  to  which  the 
model  can  handle  higher  frequencies.  For  practical  purposes 
the  PE  model  should  be  used  to  predict  propagation  loss  for 
frequencies  below  200  Hz  along  waterborne  or  shallow  bottom- 
bounce  paths  [Bef.  6]. 

The  PE  model  becomes  inappropriate  for  situations  where 
the  bottom  slopes  steeply  and  where  large  sound  speed  gradi¬ 
ents  occur  [  Bef.  6  ]. 

B.  VAIIEATIOH  AGAIIST  OOSSO  AID  CHAPMAN'S  DATA 

This  study  is  concerned  with  waterborne  paths  and 
shallow  angles.  We  are  not  specifically  interested  in 
bottom  interactions.  Even  though  the  frequencies  in  which 
we  are  interested  may  be  as  high  as  800  Hz  [Ref.  1],  we  are 
willing  to  use  computational  time  to  examine  a  full-wave 
solution.  Especially  important  for  comparison  with  Dosso 
and  Chapman's  experimental  data  is  the  fact  that  the  PE 
model  can  utilize  a  range-dependent  environment.  The  input 
parameters  described  in  Appendix  A  are  derived  from  [Ref.  1] 
and  [Bef.  7].  The  resulting  information  and  displays  show 
reasonable  agreement  with  Dosso  and  Chapman's  work. 


Naval  Ocean  Research  and  Development  Activity  (NORDA)  during 
the  ISTO’s.  A  complete  description  of  the  physics  and  algo¬ 
rithms  for  this  model  are  given  in  [Ref.  6].  The  wave  equa¬ 


tion 


[V2  +  K?  n2(r,z)]  P(r,z)  -  0 


(eqn  4.1) 


is  an  elliptic  equation  which  describes  acoustic  pressure  as 
a  function  of  range  r  and  depth  z.  It  can  be  approximated 
by  a  parabolic  partial  differential  equation  of  the  form 


f?  =  i  (  A*  B)  f» 


(eqn  4.2) 


where 


A  i  -  J- 


2>Z 


(eqn  4.3) 


B  = 


( "2  ~  I  ) 


(eqn  4.4) 


This  solution  makes  the  assumptions  that  the  sound  is 
propagating  radially  outward  from  a  source  in  a  cylindri- 
cally  symmetric  medium  and  that  the  range  r  is  very  large 
compared  to  the  size  of  the  source  and  the  wavelength  of  the 
sound:  these  are  the  cylindrical  symmetry  and  the  far- field 
assumptions.  The  source  is  further  considered  to  be  oscil¬ 
lating  harmonically  at  a  discrete  frequency.  The  ocean 
surface  is  considered  to  be  a  flat  pressure-release 
boundary,  and  the  acoustic  field  vanishes  in  the  bottom  at 
the  maximum  depth  of  the  Fourier  transform  [Ref.  6]. 

The  parabolic  equation  accounts  for  diffraction  and 
other  full-wave  effects;  it  also  allows  solution  for  range- 
dependent  environments.  The  parabolic  equation  can  be 
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Figure  4.12  Transmission  Loss  vs  Range,  630 
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Figure  4.15  Iransaission  Loss  vs  Range,  25  Hz 

(NPS  PB  Model) . 
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the  optimum  frequency  prediction  somewhat  for  the  10  June 
case. 

Putting  a-11  three  cases  together  and  looking  for  a 
trend,  we  can  tentatively  identify  a  range  of  optimum 
frequencies  around  600-900  Hz.  The  results  for  these 
shallow  sound  channels  are  not  as  clear  or  conclusive  as 
Dosso  and  Chapman's  results,  but  they  should  be  adequate  for 
the  fleet  sonar  operator  to  use  in  setting  up  search  plans. 

B.  OPTIHOB  FREQOEHCY  FOB  PROPAGATION  ACROSS  THE  LAYER 

The  second  part  of  this  investigation  involved  the  same 
sound  velocity  profiles  with  the  source  still  located  at  98 
meters  near  the  SSC  axis,  but  with  the  receiver  modeled  at 
20  meters,  well  within  the  surface  layer.  Figures  5.15 
through  5.17  show  the  results. 

Figure  5.15  for  28  June  shows  peaks  at  125  and  600  Hz  in 
the  transmission  loss  versus  frequency  curves  for  15,  24  and 
33  km.  Again  the  concept  of  cutoff  frequency  dees  not 
appear  to  have  much  meaning  in  these  cases.  Figure  5.16 
shows  peaks  at  160-200  Hz  and  at  700  Hz.  For  the  truncated 
10  June  profile.  Figure  5.17  shows  best  propagation  at 
75-100  Hz.  There  is  a  peak  at  800  Hz,  but  the  transmission 
loss  is  greater  than  at  75-100  Hz. 

Figures  5.18  through  5.20  show  relative  gain  for  our 
three  cases,  calculated  according  to  the  equation 

Relative  Gain  =TI_LAYER-  T[_ssc  lean  5.1) 

These  high- latitude  summer  profiles  do  not  show  a  cutoff 
frequency  between  the  SSC  and  the  surface  layer. 
Propagation  is  always  better  when  both  source  and  receiver 
are  contained  within  the  shallow  sound  channel  than  when  the 
source  is  in  the  SSC  and  the  receiver  is  in  the  surface 
layer. 
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Figure  5.1  High- latitude  Sound  Telocity  Profile,  10  June. 
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Figure  5.2  High- latitude  Sound  Telocity  Profile,  28  June. 
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Figure  5.4  Transaission  loss  ts  Frequency,  10  June 

(NPS  PE  Bodel)  . 
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Figure  5.19  Relative  Gain  between  Surface  Layer  and  SSC, 

10  June  (BPS  PE  Model) . 
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i. 17  Transaission  Loss  between  Surface  Layer  and  SSC, 
Truncated  Profile  (NPS  PE  Model) . 
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Trans*issi°n  Loss  between  Surface  Laver  and  SSC. 
10  June  (HPS  PE  Model).  1  ' 
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Figure  5.15  Transaission  Loss  between  Surface  Laver  and  S SC, 

28  June  (NPS  PB  Model) . 
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Figure  5.13  Transmission  Loss  vs  Range 
Truncated  Profile,  800  Hz  (NPS  PE  Model) 
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Figure  5.11  Iransaission  Loss  vs  Frequency, 
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▼I.  THE  BPS  FACT  9B  HO DEL 

A.  BOO El  DESCRIPTION 

The  Fast  Asymptotic  Coherent  Transmission  loss  mcdel, 
version  9H  (FACT  9H)  ,  like  the  split-step  Parabolic  Equation 
model,  is  accessed  through  the  Naval  Postgraduate  School* s 
IBH  3033  computer.  FACT  currently  provides  the  acoustic 
predictions  for  the  Integrated  Command  ASW  Prediction  System 
(ICAPS)  on  board  0.  S.  Navy  aircraft  carriers  and  ASW 
Operations  Centers  (ASWOCs)  worldwide.  It  is  also  the  basis 
for  a  number  of  acoustic  products  routinely  provided  by  the 
Fleet  Numerical  Oceanographic  Center,  Monterey,  to  fleet 
users.  Since  the  FACT  model  and  its  products  are  so  widely 
distributed  in  the  fleet,  it  is  of  interest  to  know  how  well 
FACT  performs  for  the  shallow  sound  channel  case. 

The  FACT  9H  model  as  installed  at  the  Naval  Postgraduate 
School  was  developed  for  the  Naval  Ocean  Research  and 
Development  Activity  (NORDA)  [Ref.  8].  It  uses  ray  tracing 
theory  with  special  corrections  to  improve  the  treatment  of 
caustics,  which  would  produce  calculations  indicating  exces¬ 
sive  acoustic  energy  at  certain  ranges  if  not  corrected. 
Unlike  the  PE  model,  FACT  9H  accepts  only  one  sound  velocity 
profile,  located  at  range  zero.  It  makes  no  allowance  for 
bottom  topography  changing  over  the  length  of  a  track.  The 
user  may  indicate  the  absorption  characteristics  of  the 
bottom,  or  specify  a  table  of  values  for  bottom  loss  versus 
grazing  angle. 

In  addition,  a  semi-empirical  expression  is  used  to 
account  for  propagation  in  a  surface  duct.  In  this  module 
the  user  may  indicate  the  wave  height,  unlike  the  NES  PE 
model  which  makes  no  provisions  for  boundary  conditions  at 
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the  surface.  As  a  ray-tracing  model,  FACT  9H  is  not 
expected  to  deal  with  waveguide  phenomena  such  as  diffrac¬ 
tion  and  leakage  as  successfully  as  the  PE  model  does.  The 
surface  duct  module  is  intended  to  remedy  this  problem  some¬ 
what  and  to  improve  the  acoustic  prediction  in  the  case  of 
surface  ducts.  No  correction  has  yet  been  applied  for 
shallow  sound  channels. 

The  investigation  of  the  FACT  9H  model  proceeded  along 
essentially  the  same  lines  as  that  of  the  PE  model.  The 
FACT  9H  model  was  first  tested  against  the  experimental 
results  cf  Dosso  and  Chapman.  Then  the  model  was  applied  to 
the  10  June  and  28  June  sound  velocity  profiles.  Input 
parameters  are  described  in  Appendix  A. 

B.  VALUATION  AGAINST  DOSSO  AND  CHAPMAN'S  DATA 

Since  FACT  9 H  does  not  allow  for  a  range-dependent  envi¬ 
ronment  or  for  wave  effects  like  diffraction  and  leakage,  it 
was  not  expected  to  perform  very  well  for  the  shallow  sound 
channel  observed  by  Eosso  and  Chapman.  Again  the  first  test 
applied  to  the  FACT  9H  output  was  a  plot  of  transmission 
loss  versus  frequency  (Figure  6.1).  As  for  the  PE  model, 
the  levels  were  obtained  by  taking  three-point  range  aver¬ 
ages  centered  at  24  km.  Since  only  the  zero  range  sound 
velocity  profile  was  used,  there  was  no  need  to  be  concerned 
about  tfce  effects  of  the  shallow  sound  channel  axis  deep¬ 
ening  over  the  eight  to  twenty  kilometer  range. 

Surprisingly,  the  transmission  loss  when  both  source  and 
receiver  were  within  the  SSC  shoved  a  gradual  decrease  and 
optimum  values  at  800-900  Hz,  consistent  with  Dosso  and 
Chapman's  actual  measurements.  Over  the  entire  frequency 
band  tested  the  FACT  9H  model  predicted  losses  between  70 
and  78  dE,  while  Dosso  and  Chapman  observed  losses  between 
about  80  and  110  dB.  The  transmission  loss  in  the  optimum 


frequency  range  is  seme  10  dB  less  than  Dosso  and  Chapman’s 
Measurements,  but  at  least  the  optimum  occurred  at  the 
proper  frequencies.  The  transmission  loss  curve  for  the 
deep  sound  channel  receiver  shoved  only  a  fairly  consistent 
increase  with  higher  frequencies. 

Figures  6.2  and  6.3  are  the  FACT  9H  transmission  loss 
curves  for  the  three  frequencies  near  optimum.  The  shallow 
receiver  curves  (Figure  6.2)  show  more  fluctuations,  prob¬ 
ably  due  to  multipath  arrivals  of  sound  reflecting  from  the 
flat  ocean  surface.  The  deep  receiver  curves  (Figure  6.3) 
are  smoothed  by  the  partial  absorption  of  sound  in  the  ocean 
bottom. 

The  plot  of  relative  gain  (Figure  6.4)  reiterates  the 
information  in  Figure  6.1  in  a  way  that  should  indicate  the 
presence  of  a  cutoff  frequency.  We  already  know  that  the 
channel  should  exhibit  a  cutoff  frequency  at  about  250  Hz, 
but  there  is  no  evidence  for  this  in  Figure  6.4.  The 
shallow  receiver  always  shows  better  propagation  than  the 
deep  receiver.  This  result  is  not  surprising  in  that  the 
FACT  9H  model  is  not  designed  to  handle  such  wave  effects. 

C.  TEE  HIGH-LATITUDE  SHALLOH  SOUHD  CHANNELS 

With  good  predictions  for  optimum  frequency  in  the  Dosso 
and  Chapman  case,  the  investigation  proceeded  to  the  10  June 
and  28  June  high- latitude,  deep  ocean  sound  velocity 
profiles.  For  the  28  June  case  (Figure  6.5)  transmission 
loss  in  the  shallow  sound  channel  did  indeed  show  a  gradual 
rise  and  fall.  The  best  frequency,  however,  was  in  the 
vicinity  of  400  Hz,  much  lower  than  the  600-900  Hz  predicted 
by  the  PE  model  (Chapter  5,  Section  A).  As  one  might 
expect,  transmission  losses  between  the  SSC  and  the  surface 
layer  are  always  significantly  greater  than  when  the  source 
and  receiver  are  both  contained  in  the  SSC. 


As  for  the  10  June  case,  FACT  9H  results  shoved  a  slight 
hut  stead;  increase  in  transmission  loss  over  the  entire 
frequency  range  (Figure  6.6) .  The  best  frequency  was,  in 
fact,  the  lowest  tested:  75  Hz.  The  transaission  loss 
curve  for  the  receiver  located  in  the  surface  layer  is  only 
a  few  dB  greater  than  the  SSC  curve.  Recall  that  the  NPS  PE 
model  also  failed  to  produce  satisfactory  results  for  the  10 
June  case.  Truncating  the  sound  velocity  profile  at  150 
meters  with  a  fully  absorbing  bottom  seemed  to  improve  the 
PE  model  prediction  somewhat,  but  Figure  6.7  does  not 
present  much  of  a  difference  for  the  FACT  9H  model. 

Perhaps  the  10  June  case  is  a  "weaker"  shallow  sound 
channel  for  modeling  purposes.  Dosso  and  Chapman's  sound 
velocity  profiles  shew  a  velocity  gradient  between  0.1436/s 
and  0.2071/s  between  the  bottom  of  the  surface  layer  and  the 
axis  of  the  shallow  sound  channel.  The  28  June  sound 
velocity  profile  has  a  velocity  gradient  of  0.1871/s.  The 
10  June  sound  velocity  profile,  however,  has  a  velocity 
gradient  of  only  0. C737/S.  There  may  be  a  critical  value 
for  the  velocity  gradient  between  the  surface  layer  and  the 
SSC  axis  that  will  produce  good  results  in  the  models. 

For  now,  however,  it  does  not  appear  that  the  FACT  9H 
model  is  successful  in  predicting  optimum  frequencies  for 
propagation  of  sound  in  shallow  sound  channels.  In  view  of 
the  extensive  usage  of  the  FACT  model  to  provide  acoustic 
transmission  loss  products  to  fleet  operations,  this  problem 
warrants  further  attention. 
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Figure  6.1  _  Transmission  Loss  vs  Frequency. 
Dosso  and  Chapman  (NPS  FACT  9H  Model). 
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Figure  6.3  Transnission  loss  between  SSC  and  DSC, 
Dosso  and  Chapean  (NPS  FACT  9H  Model) . 
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Figure  6.7  transmission  Loss  vs  Frequency. 
Truncated  Profile  ( NPS  FACT  9H  Model)  . 
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VII.  CONCLUSIONS 


A.  EE  SUITS  AID  DISCUSSION 

Two  acoustic  transmission  loss  models  were  examined  to 
determine  how  well  they  predicted  the  optimum  frequency  for 
propagation  of  sound  in  shallow  sound  channels.  The  split- 
step  Parabolic  Equation  (PE)  model  performed  satisfactorily, 
while  the  Fast  Asymptotic  Coherent  Transmission  (FACT  9H) 
model  was  less  satisfactory. 

The  EE  model  has  two  advantages:  it  allows  for  a  range- 
dependent  environment  and  it  is  a  full-wave  approximation 
for  the  wave  equation.  These  two  factors  enable  the  PE 
model  to  deal  with  wave  phenomena  like  diffraction  and 
ducting  with  a  degree  of  success.  The  PE  model  predicted 
optimum  frequencies  cf  600-900  Hz  for  the  shallow  sound 
channels  under  study,  consistent  with  the  800  Hz  actually 
measured  by  Dosso  and  Chapman.  The  model  seemed  to  produce 
tetter  results  when  the  sound  velocity  gradient  between  the 
ovetlying  surface  layer  and  the  shallow  sound  channel  was 
pronounced.  The  28  June  case,  where  the  gradient  was 
0.1371/s,  yielded  600-900  Hz  for  optimum  frequency,  but  the 
10  June  case,  where  the  gradient  was  only  0.0737/s,  did  not 
produce  an  optimum  frequency.  The  10  June  prediction  was 
improved  somewhat  by  introducing  a  fully  absorbing  bottom  at 
the  lower  boundary  of  the  shallow  sound  channel. 

The  FACT  9H  model  has  its  advantages,  too:  it  is  effi¬ 
cient  in  terms  of  computational  time  and  it  is  the  Navy's 
standard  acoustic  propagation  loss  model  at  the  time  of  this 
writing.  As  a  ray-tracing  model,  it  was  not  expected  to 
handle  situations  involving  surface  ducts  or  shallow  sound 
channels  as  successfully  as  the  PE  model,  even  with  the 
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semi-empirical  surface  duct  module  built  into  this  version 
of  FACT.  The  Dosso  and  Chapman  sound  velocity  profiles 
resulted  in  a^  surprisingly  good  prediction  of  800-900  Hz, 
which  agreed  well  with  the  experimental  results.  The  28 
June  case,  however,  predicted  an  optimum  of  about  400  Hz; 
and  the  10  June  case  did  not  predict  any  reasonable  optimum. 
Again  the  strength  of  the  sound  velocity  gradient  may 
explain  these  different  results. 

Of  our  two  choices,  the  PE  model,  despite  its  long 
computation  time,  appears  to  be  the  better  candidate  for 
predicting  optimum  frequency  of  propagation  in  a  shallow 
sound  channel.  The  FACT  model,  despite  its  widespread  use  in 
the  fleet,  is  not  a  good  predictor  for  optimum  frequency. 
The  predictions  for  tcth  models  seem  to  be  better  when  there 
is  a  large  sound  velocity  gradient  between  the  surface  layer 
and  the  axis  of  the  shallow  sound  channel. 

B.  RECGHHENDATIOHS  AID  AREAS  FOE  FURTHER  STUDY 

The  fleet  user  has  no  easy  access  to  a  Parabolic 
Equation  model  to  assess  optimum  frequencies  for  real-time 
environmental  conditions.  While  this  study  predicts  an 
optimum  range  of  600-900  Hz  where  both  source  and  receiver 
are  contained  within  80  to  110  meter  thick  channels,  it 
would  be  helpful  to  draw  up  a  table  for  other  source- 
receiver  depth  combinations  and  channels  of  varying 
thickness. 

Several  areas  merit  further  investigation.  First,  this 
study  uses  a  fully  absorbing  bottom  throughout.  This  is 
unrealistic.  It  would  be  of  interest  to  assess  bottom 
losses  more  accurately  and  determine  what  difference,  if 
any,  this  makes  in  the  prediction  of  optimum  frequency. 
Another  area  for  investigation  is  the  relative  strength  of 
shallow  sound  channels.  A  string  charnel,  where  the  sound 
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velocity  gradient  between  the  bottom  of  the  surface  layer 
and  the  axis  of  the  shallow  sound  channel,  appears  to  trap 
sound  Bore  strongly  and  to  produce  a  clearer  prediction  of 
optimuB  frequency.  It  would  be  helpful  to  find  out  how  such 
of  a  velocity  gradient  is  required  to  produce  an  optinua 
frequency  froa  each  of  the  nodels,  and  whether  the  predicted 
optiaua  varies  with  the  strength  of  the  channel. 

Finally,  other  acoustic  transmission  loss  models  are 
available.  The  implicit  finite  difference  (IFD)  PE  Bodel  is 
advertised  as  being  able  to  handle  boundary  conditions  at 
the  ocean  bottom  more  successfully  than  the  split-step  PE 
model.  The  IFD  model  should  be  investigated.  The  nost 
recent  revision  of  the  FACT  model,  version  10A,  is  nearly 
ready  for  implementation.  Technical  documents  [Bef.  9] 
indicate  that  FACT  10A  will,  among  other  things,  handle 
ducting  situations  and  boundary  conditions  better  than  FACT 
9H — this  may  affect  performance  for  shallow  sound  channels. 
A  third  type  of  acoustic  model  is  also  available — the 
PAYHCCE  model.  This  model  uses  a  combination  of  normal  mode 
and  ray- tracing  techniques  and  may  produce  better  results 
than  FACT  when  applied  to  surface  ducts  and  sound  channels. 
The  EAYMODE  model  has  recently  been  designated  as  the  Navy*s 
new  standard  acoustic  propagation  loss  model,  and  its 
performance  should  be  investigated  for  as  wide  a  variety  of 
environmental  conditicns  as  possible. 


APPEHDII  1 

II PUT  PARAMETERS  TOR  THE  IPS  PE  AID  FACT  9H  MODELS 


Table  I  describes  the  input  options  available  to  the 
user  through  the  Naval  Postgraduate  School*s  split-step 
Parabolic  Equation  (PE)  and  Fast  Asymptotic  Coherent 
Transmission  loss  (FACT  9H)  models.  Tables  II  through  VI 
are  the  sound  velocity  profiles  used  in  this  study. 


TABLE  I 

Inpat  Parameters  for  the  IPS  PE  and  FACT  9H  Models 


PAHAMITEES 

PE  MODEL 

FACT  9  H 

Input  units 

English 
or  metric 

English 
or  metric 

Source  depths 

One 

One 

Beceiver  depths 

Maximus  20 

One 

Frequencies 

One 

Maximum  6 

Maxiaun  range 

Unlimited 

999  kilometers 

Bange  step 

Minimum  .01  nm 
or  selected  by 
program 

Selected  by 
program 

Humber  of  SVPs 

Maximum  50 

One 

pepth  levels 
in  SVPs 

Maximum  50 

Maximum  50 

Bottom  less 

1)  Fully  absorbing 

1)  Bottom  type 

0  (reflecting) 
to  9  (absorbing) 

2)  User’s  loss  vs 
grazing  angle  table 

2)  Usep’s  loss 
vs  grazing  angle 
table 

31  User’s  depth  vs 
attenuation  table 

3)  Internal  loss 
vs  grazing  angle 
table 

Variable 
bottoa  depth 

Maximum  100  ranges/ 
depths 

No 

Correction  for 
spherical  earth 

Available 

Available 

Correction  for 
for  volume 
attenuation 

Available 

Available 

Vertical  beam 
size  of  source 

1-33  degrees 

No 

Special  handling 
for  surface  duct 

Not  required 

Yes 

Have  height 

No 

Yes 

Critical  angle 

Not  reguired 

1)  Preset  value 
28.65  degrees 

2)  User  input 


TIBLE  II 

Dosso  and  Chap nan's  Sound  Telocity  Profile  at  0  Kiloneters 


DEPIH(FT) 

C (FT/S) 

DEPTH  (M) 

C  (M/S) 

0. 

4927.82 

0. 

1502.00 

45. 

4927.82 

14. 

1502.00 

66. 

4922.90 

20. 

1500.50 

87. 

4906.50 

26. 

1495.50 

115. 

4890.75 

35. 

1490.70 

137. 

169. 

42. 

52. 

U»:38 

198. 

4873.36 

60. 

1485.40 

232. 

4867.45 

71. 

1483.60 

251. 

4866.47 

76. 

1483.30 

273. 

4865.48 

83. 

1483.00 

291. 

4864.83 

89. 

1482.80 

310. 

4864.83 

95. 

1482.80 

325. 

4864.50 

99. 

1482.70 

338. 

4864.50 

103. 

1482.70 

357. 

4863.52 

109. 

1482.40 

369. 

4863.84 

112. 

1482.50 

384. 

4864.83 

117. 

1482.80 

403. 

4865.48 

123. 

1483.00 

433. 

4866.14 

132. 

1483.20 

476. 

510. 

4867.12 

4867.12 

145. 

155. 

1483.50 

1483.50 

543. 

4866.47 

166. 

1483.30 

589. 

4866.47 

179. 

1483.30 

690. 

4862.20 

210. 

1482.00 

779. 

856. 

4859.58 

4858.59 

238. 

261. 

1481.20 

1480.90 

946. 

4857.61 

288. 

1480.60 

1026. 

4854.66 

313. 

1479.70 

1115. 

4854.00 

340. 

1479.50 

1312. 

4849.41 

400. 

1478. 10 

1640. 

4849.08 

500. 

1478.00 

1969. 

4849.73 

600. 

1478.20 

2625. 

4852.36 

800. 

1479.00 

3281. 

48  56.62 

1000. 

1480.30 

3937. 

4860.89 

1200. 

1481.60 

4921. 

4870.08 

1500. 

1484.40 

6562. 

4891.07 

2000. 

1490.80 

8202. 

4916.01 

2500. 

1498.40 

See  Figure  3.1. 


TABLE  III 

Dosso  and  Chapnan* s  Sotmd  Telocity  Profile  at  16  Kilonetexs 


DEPTH  (FT) 


C{  FT/S) 


DEPTH  (M) 


C  (H/S) 


m. 

1148. 

1181. 

1214. 

1247. 

1280. 

1312. 

1345. 

1411. 

m. 

1640. 

1969. 

2625. 

1114: 

4921. 

6562. 

8202. 

9514. 


4923.  55 

4923.88 

4926.51 
4924.21 

4924.21 

4898.62 
4884.18 

4879.26 

4877.62 

4873.69 

4867.12 

4864.50 

4865.48 

4863.84 

4865.16 

4863.52 

4860.89 
4860.56 

4861.22 
4862.86 

4861.55 

4861.55 

4860.23 
4859.91 
4859.91 

4858.59 

4858.27 
4857.61 
4855.64 
4853.67 

4852.69 

4850.72 
4851.37 
4850.39 
4850.39 

4349.73 
4849.41 
4850.39 
485C.72 
4852.69 

4870.08 

4891.07 

4916.01 

4937.01 


80. 

90. 

100. 

110. 

120. 

130. 

140. 

150. 

160. 

180. 

210. 

240. 

260. 

280. 

300. 


1  oog. 
1200. 
1500. 
2000. 
2500. 
2900. 


1500.70 

1500.80 

1501.60 
15C0.90 

1500.90 

1811:18 

1487.20 

1486.70 

1485.50 

1483.50 

1482.70 
1483.00 

1482.50 

1482.90 

1482.40 

1481.60 

1481 .50 

1481.70 

1482.20 

1481.80 

1481.80 

1481.40 
1481.30 

1481.30 

UfS:f8 

1480.60 
1480.00 

1479.40 

1479. 10 

1478.50 

1478.70 

1478.40 
1478.40 

1478.20 

1478.  10 

1478.40 
1478.50 

1479.  10 

1480.30 

1481 .60 

1484.40 

1490.80 

1498.40 

1504.80 


See  Figure  3.2. 
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Dosso  and  Chap «an's  Sound  Telocity  Profile  at  20 
and  33  Kiloneters 


DEPTH  (FT) 

C ( FT/S) 

DEPTH  (ft) 

C  (H/S) 

0. 

4922.57 

0. 

1500.40 

33. 

4921.26 

10. 

1500.00 

66. 

48  96.32 

20. 

1492.40 

98. 

4885.82 

30. 

1489.20 

131. 

4878.28 

40. 

1486.90 

164. 

4875.66 

50. 

1486.10 

197. 

4875.33 

60. 

1486.00 

230. 

4873.69 

70. 

1485.50 

262. 

4871.06 

80. 

1484.70 

279. 

4866.47 

85. 

1483.30 

295. 

4864.50 

90. 

1482.70 

312. 

4863.52 

95. 

1482.40 

328. 

4864.17 

100. 

1482.60 

344. 

4863.84 

105. 

1482.50 

361. 

4863.84 

110. 

1482.50 

394. 

4865. 16 

120. 

1482.90 

427. 

4864.83 

130. 

1482.80 

459. 

4865.81 

140. 

1483.10 

492. 

4865.48 

150. 

1483.00 

558. 

4865.81 

170. 

1483.10 

591. 

4864.17 

180. 

1482.60 

623. 

4863.52 

190. 

1482.40 

656. 

4860.56 

200. 

1481.50 

722. 

4860.23 

220. 

1481.40 

787. 

4857.61 

240. 

1480.60 

853. 

4856.30 

260. 

1480.20 

919. 

4853.67 

280. 

1479.40 

984. 

4850.72 

300. 

1478.50 

1083. 

4849.73 

330. 

1478.20 

1148. 

4851.05 

350. 

1478.60 

1312. 

4848.75 

400. 

1477.90 

1640. 

4849.08 

500. 

1478.00 

1969. 

4849.41 

600. 

1478.10 

2625. 

4851.70 

800. 

1478.80 

3281. 

4855.97 

1000. 

1480. 10 

3937. 

4860.89 

1200. 

1481.60 

4921. 

4870. 08 

1500. 

1484.40 

6562. 

4891.07 

2000. 

1490.80 

8202. 

49  16.01 

2500. 

1498.40 

8858. 

4927.82 

2700. 

15C2.00 

TABU  T 

High- latitude  Sound  Telocity  Profile  for  10  June 


DEPTH  (FT) 

C (IT/S) 

DEPTH 

0. 

47S6.59 

0. 

131. 

47S6.59 

40. 

164. 

4780.  18 

50. 

328. 

4780. 18 

100. 

443. 

4773.62 

135. 

ill: 

4798.23 

484*.  52 

1640. 

4824.47 

500. 

2281. 

4847.11 

1000. 

6562. 

4890.42 

2000. 

S843. 

4943.57 

3000. 

3123. 

4999.34 

4000. 

See  Figure  5. 1 . 


C  (M/S) 


1462.00 
1462.00 
1457.00 
1457. 0C 
1455.00 


14 

1470.50 

1477.40 

1490.60 


1506.80 

1523.80 


TABLE  VI 

High- latitude  Sound  Telocity  Profile  for  28  June 


EEFTH  (FT) 

C ( FT/S) 

0. 

48C6.43 

115. 

4804.46 

164. 

4775.26 

262. 

4771.65 

328. 

4765.42 

459. 

4799.54 

656. 

4811.68 

1640. 

4  827.43 

3281. 

4847.11 

6562. 

4890.42 

iSHS: 

!!&!Z 

DEPTH  (N) 

C(H/S) 

0. 

1465.00 

35. 

1464.40 

50. 

1455.50 

80. 

1454. 4C 

100. 

1452.50 

140. 

1462.90 

200. 

1466.60 

500. 

1471.40 

1000. 

1477.40 

2000. 

1490.60 

3000. 

4000. 

1121:18 

See  Figure  5.2 
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